Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1582 1652 (N46W53 L=256)  1659    226        3 (1481/1606-1687/1693)  -6.1 W41-W70 N41-N49
      2 1023 1562 (N38E02 L=201)   456    105        1 ( 972/1539-1071/1588)  18.4 E05-W01 N35-N39
      3  325 1505 (N30E62 L=141)   766     85        1 ( 289/1475- 374/1525)   4.4 E68-E55 N27-N31
      4 1090 1482 (N32W03 L=206)   360     69        1 (1064/1438-1114/1517) -59.0 W01-W04 N29-N34
      5  496 1492 (N30E42 L=160)   311     58        0 ( 460/1476- 524/1508)  11.3 E45-E40 N29-N31
      6 1081 1345 (N23W03 L=206)   357     73        0 (1057/1307-1108/1380)  68.2 W01-W04 N20-N25
      7  179 1166 (N08E66 L=137)   719    108        0 ( 161/1111- 198/1218) -40.6 E70-E62 N04-N10
      8 1490 1089 (N08W30 L=232)  1364    295        0 (1402/ 958-1543/1202) -60.3 W23-W33 N00-N14
      9 1637 1065 (N06W41 L=243)   421    103        2 (1609/1034-1674/1115) -59.0 W38-W43 N04-N09
     10 1720  813 (S09W49 L=252)  1142    156        2 (1677/ 742-1757/ 885)  66.0 W46-W52 S13-S04
     11  196  842 (S12E63 L=140)   432     38        1 ( 176/ 821- 219/ 862) -33.7 E65-E59 S13-S11
     12 1104  761 (S13W06 L=208)   200     50        0 (1082/ 736-1130/ 789) -45.0 W04-W07 S14-S11
     13  187  731 (S20E68 L=134)   874     84        1 ( 162/ 694- 206/ 772)  31.0 E76-E64 S22-S16
     14 1634  600 (S23W46 L=249)  1383    170        3 (1574/ 545-1690/ 663)  53.1 W41-W50 S26-S18
     15  685  607 (S25E22 L=181)  1151    204        2 ( 612/ 566- 776/ 660) -10.3 E26-E15 S27-S21
     16  493  598 (S26E37 L=165)   285     80        0 ( 454/ 580- 520/ 631)  33.7 E40-E34 S27-S23
     17 1161  217 (S56W19 L=222)  2070    186        3 (1074/ 189-1257/ 257)   8.5 W09-W29 S59-S50
  

    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