Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  207 1402 (N28E73 L=328)  1085    102        1 ( 185/1366- 235/1439) -51.3 E80-E68 N25-N30
      2 1592 1313 (N19W41 L= 82)   370     86        0 (1557/1280-1626/1345) -51.3 W38-W43 N16-N21
      3  436 1289 (N21E41 L=  0)   405     79        1 ( 400/1262- 464/1319)  31.0 E43-E38 N19-N23
      4  318 1284 (N21E52 L=349)   513     65        0 ( 282/1256- 355/1313) -18.4 E56-E47 N19-N22
      5 1830 1272 (N15W64 L=105)   984     73        1 (1801/1232-1855/1311) -68.2 W60-W67 N12-N17
      6  664 1162 (N13E23 L= 19)   246     70        0 ( 634/1147- 696/1178) -14.0 E24-E20 N12-N14
      7 1374 1041 (N04W22 L= 63)   447    168        0 (1315/ 984-1441/1092)  33.7 W17-W26 N00-N06
      8 1692 1015 (N00W45 L= 86)   435    165        0 (1681/ 931-1705/1092)  83.7 W44-W46 S04-N05
      9  476  972 (N02E36 L=  6)   487    221        0 ( 402/ 911- 562/1025) -22.6 E41-E29 S01-N05
     10 1369  812 (S10W21 L= 62)  1353    534        0 (1188/ 683-1481/ 995)  45.0 W09-W28 S17-N00
     11 1036  920 (S02W00 L= 42)   263     91        0 (1004/ 882-1064/ 966)  53.1 E01-W02 S04-N00
     12 1846  928 (S06W61 L=102)   502     56        0 (1809/ 914-1876/ 944)  -7.1 W56-W64 S07-S05
     13  664  842 (S06E23 L= 18)   512    178        0 ( 598/ 812- 711/ 906)  35.5 E27-E20 S07-S02
     14  308  773 (S11E52 L=349)   310     44        0 ( 283/ 745- 328/ 796) -45.0 E54-E50 S12-S09
     15  432  654 (S18E43 L=358)   348     89        1 ( 405/ 616- 473/ 680)  33.7 E45-E38 S20-S16
     16  360  651 (S18E50 L=351)   354     58        1 ( 339/ 620- 384/ 677) -63.4 E51-E46 S20-S16
  

    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