Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1431 1517 (N35W32 L= 89)   457     68        0 (1406/1479-1453/1559) -90.0 W30-W33 N31-N37
      2 1657 1490 (N31W51 L=108)   790    114        1 (1621/1438-1704/1532) -60.3 W47-W53 N27-N34
      3  311 1480 (N32E58 L=359)   829    218        1 ( 277/1456- 366/1513) -15.3 E62-E51 N30-N34
      4  165 1224 (N15E69 L=348)   581     84        1 ( 148/1202- 185/1247) -21.8 E71-E65 N13-N15
      5  734 1181 (N14E18 L= 39)   217     62        1 ( 700/1153- 764/1208) -45.0 E20-E16 N11-N15
      6 1133  970 (N00W07 L= 64)   320     98        0 (1084/ 950-1188/ 991)   8.1 W03-W10 N00-N01
      7 1857  887 (S08W63 L=121)  1479    284        2 (1826/ 789-1915/ 977)  45.0 W59-W72 S13-S02
      8  539  753 (S13E32 L= 25)   410    103        0 ( 481/ 726- 580/ 783) -23.2 E36-E29 S14-S11
      9  359  615 (S23E50 L=  7)  1531    377        1 ( 286/ 536- 424/ 704)  35.8 E60-E42 S29-S16
     10 1041  573 (S25W01 L= 58)   330     64        1 (1000/ 560-1078/ 588) -11.3 E02-W03 S25-S24
  

    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