Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  972 1871 (N46E37 L=  8)  1017    111        3 ( 933/1839-1013/1901)  80.5 E40-E33 N41-N49
      2 1075 1698 (N35E19 L= 27)  2018    472        7 (1013/1545-1121/1854) -42.4 E32-E09 N24-N46
      3  371 1585 (N10E62 L=343)  2192    485        5 ( 318/1458- 420/1698) -67.6 E69-E53 N01-N18
      4  700 1399 (N06E29 L= 16)   867    236        1 ( 628/1340- 792/1479) -11.3 E36-E21 N03-N09
      5  930 1316 (N07E13 L= 32)   229     52        1 ( 911/1290- 955/1342)  71.6 E14-E11 N05-N08
      6 1037 1317 (N10E07 L= 38)   251     74        0 ( 994/1301-1077/1331)   0.0 E09-E04 N08-N10
      7  506 1311 (S04E39 L=  7)   329    115        0 ( 457/1283- 546/1329)   0.0 E42-E34 S04-S02
      8 1097 1243 (N07E02 L= 43)   206     32        0 (1074/1221-1116/1265)  71.6 E02-W00 N05-N08
      9 1006 1095 (S04E03 L= 42)   803    258        2 ( 882/ 993-1119/1179)  61.9 E07-W01 S12-N03
     10  560 1003 (S20E27 L= 18)   207     41        0 ( 529/ 985- 589/1022)   0.0 E29-E24 S20-S18
     11 1868  925 (N14W59 L=104)   552    108        1 (1854/ 889-1895/ 973) -59.0 W56-W63 N11-N16
     12  529  914 (S26E27 L= 18)   325     53        0 ( 507/ 884- 552/ 941)  90.0 E29-E25 S27-S23
  

    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