Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1599 1531 (N33W48 L= 40)   537    122        0 (1553/1504-1650/1559) -33.7 W44-W52 N31-N35
      2  749 1344 (N25E18 L=334)  3720   1067        7 ( 338/1105-1072/1548)  -7.9 E57-W03 N09-N39
      3 1481 1421 (N27W34 L= 26)   301     79        1 (1444/1407-1521/1440) -31.0 W32-W37 N25-N28
      4 1693 1235 (N13W48 L= 40)   990    251        2 (1612/1146-1783/1318) -42.7 W42-W55 N07-N19
      5  347 1170 (N14E47 L=304)   651    169        1 ( 328/1095- 372/1259) -78.7 E49-E44 N09-N19
      6  168  924 (S01E67 L=284)   814    182        1 ( 130/ 904- 214/ 942)   8.1 E74-E60 S03-N00
      7 1134  908 (S03W06 L=358)   196     38        1 (1109/ 890-1164/ 923)   0.0 W04-W08 S04-S02
      8 1425  861 (S08W25 L= 16)   565    141        1 (1362/ 817-1482/ 913) -32.0 W20-W28 S10-S03
      9  928  688 (S16E07 L=344)  1015    328        1 ( 805/ 576-1027/ 815)  43.2 E16-W00 S23-S08
     10 1574  761 (S15W36 L= 28)   377     76        1 (1537/ 743-1608/ 783)  11.3 W33-W38 S15-S13
     11 1134  522 (S29W06 L=357)  1964    481        3 ( 974/ 384-1326/ 689)  41.0 E06-W18 S38-S17
  

    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