Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1400 1560 (N32W38 L= 33)   441    168        1 (1321/1539-1464/1586)  -9.5 W31-W43 N30-N33
      2 1086 1429 (N30W11 L=  6)   719    153        1 (1020/1379-1170/1467)   9.5 W05-W17 N27-N31
      3 1659 1352 (N14W49 L= 44)  1266    202        1 (1587/1287-1733/1394) -45.0 W44-W54 N08-N18
      4 1126  961 (N01W05 L=360)  1781    461        8 (1027/ 771-1175/1171) -85.4 W01-W07 S10-N15
      5  515 1068 (N16E32 L=323)   515    122        0 ( 464/1042- 568/1105) -23.2 E35-E28 N13-N18
      6 1431  901 (S08W23 L= 18)   798    196        2 (1353/ 841-1520/ 950) -42.0 W18-W28 S12-S03
      7 1302  799 (S12W13 L=  8)   387    104        1 (1283/ 745-1322/ 860)  63.4 W11-W15 S14-S08
      8 1117  734 (S13W01 L=356)   217     55        1 (1088/ 708-1143/ 759)  45.0 W00-W02 S13-S11
      9 1308  581 (S26W11 L=  6)   650    122        1 (1277/ 522-1340/ 640)  54.5 W08-W14 S29-S22
     10  962  350 (S37E17 L=338)  1247    308        3 ( 884/ 278-1011/ 466)  40.9 E26-E11 S41-S28
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y